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Abstract 

Searches for the permanent electric dipole moments (EDMs) of diamagnetic atoms provide pow- 
erful probes of CP-violating hadronic and semileptonic interactions. The theoretical interpretation 
of such experiments, however, requires careful implementation of a well-known theorem by Schiff 
that implies a vanishing net EDM for an atom built entirely from point-like, nonrelativistic con- 
stituents that interact only electrostatically. Any experimental observation of a nonzero atomic 
EDM would result from corrections to the point-like, nonrelativistic, electrostatic assumption. We 
reformulate Schiff 's theorem at the operator level and delineate the electronic and nuclear operators 
whose atomic matrix elements generate corrections to "Schiff screening" . We obtain a form for the 
operator responsible for the leading correction associated with finite nuclear size - the so-called 
"Schiff moment" operator - and observe that it differs from the corresponding operator used in 
previous Schiff moment computations. We show that the more general Schiff moment operator 
reduces to the previously employed operator only under certain approximations that are not gen- 
erally justified. We also identify other corrections to Schiff screening that may not be included 
properly in previous theoretical treatments. We discuss practical considerations for obtaining a 
complete computation of corrections to Schiff screening in atomic EDM calculations. 
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I. INTRODUCTION 



The search for CP-violation (CPV) in and beyond the Standard Model (SM) is important 
to particle, nuclear, and atomic physics and their intersections with cosmology. The CPV 
in the electroweak sector of the SM, parametrized by the complex phase in the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix, adequately accounts for experimental observations of 
CPV effects in neutral kaon and B meson systems. The results of searches for the permanent 
electric dipole moments (EDMs) of neutron (dn) and diamagnetic atoms (0?^) imply that 
CPV in the strong sector of the SM, parametrized by the so-called "^-term" are vanishingly 
small: |^| < 10~^° [1, 2]. The so-called strong CP-problcm associated with the unnaturally 
small value of 6 may find resolution in spontaneously broken Peccei-Quinn symmetry if the 
corresponding Goldstone boson - the axion - is found. If so, strong CPV may hold the key 
to one of the outstanding puzzles in cosmology, as the axion is a viable candidate for particle 
dark matter. 

Nevertheless, another cosmological problem - the origin of the non-vanishing baryon 
asymmetry of the universe (BAU) - would remain unsolved. As noted by Sakharov nearly 
four decades ago [3] , if the initial conditions of the universe were matter-antimatter sym- 
metric, the particle physics of the evolving universe would have to contain CPV interactions 
(along with baryon number violation and a departure from thermal equilibrium) in order to 
produce a BAU. It has been subsequently noted that the CPV contained in the SM is not 
sufficiently effective to produce the observed BAU, so an explanation requires new sources 
of CPV. If such interactions involve particles having masses < a few TeV, then current 
searches for atomic, hadronic, and leptonic EDMs could soon yield nonzero results. Exist- 
ing limits already provide stringent constraints on new electroweak-scale CPV, and the next 
generation of experiments are poised to push the sensitivity to new CPV sources by several 
orders of magnitude. If new electroweak-scale CPV in models like supersymmetry is, indeed, 
responsible for the BAU, then one would expect nonzero results in this new round of EDM 
searches (see, e.g., Ref. [4] and references therein. For a recent review of EDMs and CPV, 
see Ref. [5]). 

If a nonzero EDM is observed, then it will take a set of complementary EDM experiments 
- together with the results of future collider searches and precision electroweak studies - to 
identify the nature of the CPV. In this respect, searches for EDMs of diamagnetic atoms 
play a key role. Experimentally, new efforts to probe with Ra [6, 7], Xe[8], and Rn [9] 
are presently underway. The sensitivity of these atomic EDMs to either |^| or CPV beyond 
that of the SM could exceed that obtained for the ^^^Hg atom [2] , which currently sets the 
standard for dA, by several orders of magnitude. The ^^^Hg standard is also expected to 
improve by an order of magnitude, when current experiments are completed [10]. Theoreti- 
cally the interpretation of dA is challenging because of the interplay of atomic, nuclear, and 
hadronic physics. Specifically, it has long been recognized that an EDM residing on the nu- 
cleus of a neutral diamagnetic atom would not be detectable, apart from small corrections: 
Schiff [11] showed that any neutral system of electrically charged, nonrelativistic, point-like 
constituents interacting only electrostatically will have no net EDM. ^ 

Schiff's result can be understood at the classical level. An atomic EDM is probed by 
placing the atom in a combination of external magnetic and electric fields. If the atom has 



In fact, part of this shielding effect has already been pointed out by Purcell and Ramsey [12] much earlier, 
but it is Schiff who first clearly demonstrates this effect in quantum mechanics. 
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a net EDM, then its interaction with the external electric field Eg^t would lead to a shift in 
the Larmor precession frequency that depends linearly on Ecxt- However, a neutral system 
in the presence of an external electric field cannot accelerate. Thus, the atomic constituents 
- assumed here to be point-like - must rearrange so as to screen out the effect of E'ext 
at the location of each of the charged constituents. Consequently, if any of the charged 
constituents (electrons or nucleus) possesses a nonzero EDM, changes in the interaction 
energy due to E^^t will be exactly cancelled by those from the internal fields induced by 
atomic polarization. To the extent that the interaction of the atom with E^xt arises solely 
from the EDMs of its point-like, nonrelativistic constituents, there will be no shift in the 
Larmor frequency linear in Eext and, hence, no effective atomic EDM. Corrections to this 
exact screening will arise if any of the constituents having an EDM (a) moves relativistically, 
leading to a breakdown of the classical picture, (b) has a finite size, leading to a breakdown 
of the point-like assumption, or (c) has non-electrostatic interactions that become modified 
in the presence of the atomic rearrangement needed to achieve electrostatic equilibrium. In 
addition, CPV interactions between constituents can lead to an atomic EDM that evades 
the Schiff screening. 

Because of such corrections and the exquisite sensitivity of the measurements, atomic 
EDM experiments place important constraints on the electron EDM and on CPV interactions 
in the nucleus. The present limits on the electron EDM d^, for example, are derived from 
heavy paramagnetic atoms, where the Schiff theorem is evaded by the rclativistic motion of 
the atomic electrons. In the case of diamagnetic atoms, nonzero effects due to hadronic CPV 
arise because of corrections (b) and (c) , specifically the finite size of the nucleus that leads 
to imperfect shielding and magnetic and other, higher-order multipole interactions between 
the nucleus and atomic electrons. The finite-size correction becomes much more important 
in heavier atoms: to the extent that the atomic electrons penetrate the atomic nucleus, a 
residual CPV interaction arises. This penetration produces an atomic EDM proportional to 
the so-called "Schiff moment" (S) of the nucleus, often expressed as an r^-weighted integral 
over the difference in the nuclear charge and EDM distributions. Additional contributions 
to the atomic EDM arc generated by hypcrfinc (non-electrostatic) interactions between 
atomic electrons and CP-conserving higher-order nuclear moments, such as the Mi (magnetic 
dipole) and C2 (charge quadrupole), as well as higher-order CPV nuclear moments that 
couple to spatially varying fields, such as the M2 (magnetic quadrupole). 

These general features of Schiff screening and their corrections were identified in Schiff's 
original paper, and subsequently discussed in other work by Sandars [13], Feinberg [14], and 
Sushkov, Flambaum, and Khriplovich [15] (for later discussions, see, e.g., Refs. [16-20]). The 
formulation of the corrections to Schiff screening in Refs. [14, 15] concentrated on the Schiff 
moment effect, and approximate expressions for S were derived. The latter have been used 
in subsequent theoretical treatments of the EDMs of ^^^Hg and other atoms, from which 
limits on CPV parameters, such as \9\, have been derived (see, e.g. Ref. [21] and references 
therein). While these treatments are correct at a qualitative level, more precise derivations 
are possible, and likely important given the prospect that nonzero atomic EDMs will soon be 
seen. In what follows, we provide a systematic derivation of Schiff's theorem and associated 
corrections in a spherical multipole formalism. When comparisons are made to earlier work, 
several refinements are apparent: 

(i) Our approach expresses the electron-nucleus interaction in terms of spherical multipole 
operators, separating contributions associated with electron penetration from those 
associated with a point-like nucleus. In so doing, we obtain a form for the atomic 
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operator O^^^^^ that describes the leading corrections to Schiff screening associated 
with electronic penetration into the nucleus. This atomic operator can be applied to 
any atomic EDM and gives a more general characterization of the leading finite size 
corrections than has been used in previous studies. 

(ii) We subsequently derive the nuclear Schifi' moment operator S by considering matrix 
elements of OgJ;^^^"^ between electronic S- and P-states as is relevant for many atomic 
EDMs. Our result for S, the analog of the Schiff moment operator used in many 
earlier studies, has a more general form than the one in common use. We discuss the 
approximations under which our expression for S reduces to the form conventionally 
employed, and observe that these approximations are not generally justified. We 
advocate the use of the operator S derived below that does not require adoption 
of any approximations. We also note that the operator S derived for electronic S- 
P transitions is not appropriate for other transitions, such as the close-lying D-P 
transition in a metastable Ra atom where a large enhancement factor is found [22] . 

(iii) We derive expressions for the operators that characterize interactions between the po- 
larized atomic cloud and CP-conserving nuclear moments, such as the magnetic dipole 
and charge quadrupole moments. We also include finite-size, electronic penetration 
effects associated with these multipole moments. 

(iv) Wc obtain additional electronic penetration corrections involving higher CPV multi- 
pole moments of the nuclear charge and current densities, such as the charge octupole 
moment and magnetic quadrupole moment. 

(v) We revisit corrections to Schiff screening of electron EDMs and include corrections 
associated with both the external and nuclear vector potentials that have not been 
included in earlier studies. 

In addition to deriving the operators associated with the foregoing effects, we also discuss 
practical considerations involved in computing atomic EDMs. For example, a complete 
computation of effects associated with hadronic CPV would require simultaneous diagonal- 
ization of the electronic and nuclear Hamiltonian, taking into account effects such as the 
nuclear response to the polarized atomic cloud and the corresponding interaction with the 
external electric field. Practically speaking, carrying out such diagonalization of the full 
space of atomic wave functions is not possible even in the limit of purely CP-conserving 
interactions. In typical atomic calculations, the nuclear charge and current distributions are 
taken as c-number sources to which the atomic electrons respond, and "back reactions" of 
the nuclear charge and current distributions to the atomic electrons - the so-called nuclear 
polarization corrections - are neglected. This procedure amounts to treating the atomic 
states as direct products of the nuclear ground state and different electronic states, with the 
properties of the nuclear ground state acting as a c-number input for diagonalization of the 
electronic Hamiltonian: 

l^)atom = lg-S-)Ar® l^)e > (l) 

where the A/" and e subscripts refer to the nuclear ground state and electronic states (labeled 
by quantum numbers k), respectively. 
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One would naively expect the nuclear polarization corrections to atomic properties asso- 
ciated with this factorization to be small, going as powers of the nuclear radius, Rj^f., divided 
by the effective Bohr radius, ao- For electronic S-states, however, the situation is more sub- 
tle. Theoretical studies of muonic atoms, for which the effective Bohr radius is ~ 200 times 
smaller than for ordinary atoms, indicate that corrections to the energy splitting between 
the 2Si/2 and 2P3/2 states depends on the square of the muon wave function at the origin 
and an integral over the nuclear photoabsorption cross section [23-25]. Numerically, the 
theoretical polarization correction to A£'(2P3/2 — 2S1/2) in (/x^He)"*" is about 0.2% of the 
measured splitting. In the case of (e^He)+, the authors of Refs. [23, 24] find that the nuclear 
polarization correction to the S-P splitting is a few xlO~^ times smaller in magnitude than 
the corresponding finite size correction, which depends on the square of the electronic wave 
function at the origin times the mean square nuclear charge radius. 

The Schiff moment interaction itself characterizes finite-size corrections of order {Rj^/aoY 
relative to the energy associated with the interaction of the atomic electrons with the EDM 
of a point-like nucleus. Based on the studies of Refs. [23-25], one would expect the nuclear 
polarization corrections to the finite-size effects characterized by the Schiff moment to be 
negligible. In the present analysis, however, we have not quantified these polarization effects, 
and leave this task to future work. Instead, we will proceed by formulating the Schiff 
theorem at the operator level as far as possible without making explicit reference to the 
atomic states or invoking the factorization approximation, thereby avoiding the issue of 
nuclear polarization corrections. Only when attempting to compare our Schiff moment 
operator with the corresponding operator used previously in the literature will we adopt the 
factorization ansatz of Eq. (1). Even in this case, we find substantial differences with Schiff 
moment operator used in earlier analyses. 

With the foregoing caveats in mind, we believe that our reformulation of Schiff 's theorem 
represents a useful refinement that will allow for a more complete inclusion of hadronic 
and nuclear structure contributions to dA and provide for a sharper confrontation between 
theory and experiment. We organize our reformulation in the remainder of the paper as 
follows. In Section 11 we give our multipolc operator formulation of the theorem, including 
a delineation of the various operators that characterize corrections to Schiff screening. For 
completeness, we include here a discussion of both the hadronic CPV effects as well as 
those associated with electronic EDMs. We include here for the first time the corrections 
associated with non-electrostatic interactions between atomic electrons and both nuclear and 
external sources. In Section III, we obtain the full atomic Schiff moment operator Ogchig'^ 
that is independent of the atomic states and use it, together with the factorization ansatz, 
to derive an effective nuclear Schiff moment operator S that can be compared with the 
operator used previously in the literature. We also derive the operator that characterizes 
the leading corrections to Schiff screening associated with magnetic interactions between 
the electrons and nucleus. In Section IV, we summarize our results and discuss the practical 
implementation of Schiff's theorem including the issues alluded above. A few technical 
details are given in the Appendices. 



II. SCHIFF'S THEOREM: A MULTIPOLE OPERATOR FORMULATION 

The interaction of an external field with the atomic EDM is odd under both parity (P) 
and time- reversal (T) transformations but even under charge conjugation (C). Assuming 
CPT invariance, as we do throughout this paper, the observation of a nonzero d^-induced 
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frequency shift is, thus, equivalent to the observation of CP-violation. Henceforth, we will 
refer to atomic, electronic, and nuclear moments that violate both P and T symmetries as 
PVTV moments, while those that respect these symmetries are PCTC moments. Effects 
that violate P but conserve T, such as the nuclear anapole moment, are also of considerable 
interest, but we do not consider these effects here (for recent reviews, sec Rcfs. [26. 27]). 
Effects that conserve P but violate T, which in combination with the weak interaction can 
induce PVTV interactions, are also not explicitly treated here. 



A. The atomic Hamiltonian 

To evaluate the consequences of PVTV interaction of atomic properties, we separate the 
complete atomic Hamiltonian into PCTC (Hq) and PVTV (Hj) terms, treating the latter 
as a perturbation: 

-f^atom = Hq + Hj . (2) 

We also include the interaction of an external field with the atom in Hj. ^ Before identifying 
the specific operators in Hq and Hj, it is useful to delineate the different interactions between 
atomic electrons, the electrons and nucleus, and external fields and the atomic constituents: 

1. The interactions between atomic electrons is 

f7-(ee) _ T/(ee) , fr{ee) /o\ 

where we adopt a notation (used throughout this paper) where the superscript (ee) 
denotes a PCTC electron-electron interaction with PCTC couplings on both electrons, 
while (ee) denotes a PVTV interaction that arises from a PCTC coupling on one 
electron and a PVTV coupling on the other. Specifically, V^^^^ denotes the PCTC 
Coulomb and Breit interactions between electrons 

Kif = ? E (4" - ■ ^:'") (4) 



and V^nt'^'' denotes the PVTV electron-electron interaction, which we write as 

V^'^ = f E 4 /? [cr. • El^^ + z a. • )] + . . . . (5) 

1=1 

The term given explicitly is the interaction of each electron's EDM with the electric 
and magnetic fields (^^I^^S^^-*) created by the other electrons. In addition, the " 
indicates that possible additional terms could be added to this, new exchanges gener- 
ating PVTV interactions. Here, a is the fine structure constant and Pi and cxi are the 
Dirac matrices acting on the ith electron. The scalar potential = (j)^^\xi) is the 
potential exerted by the other Z — 1 electrons at the position Xi of the ith electron; 
the notation is similar for the vector potential A'f\ and for the resulting electric and 
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As we are interested on permanent EDMs instead of induced EDMs, it is the first-order Stark effect that 
one is after. Therefore, the weak external field limit is adequate. 
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magnetic fields E^;^^ and B^^K The 1/2 factor in the electron-electron interaction is 
introduced to ensure that one sums over distinct pairwise interactions. ^ 

2. The interaction between the electrons and nucleus (jV) is 



-«^[r>-a..Af)]+v;2'^, (6) 



— — ut 

i=l 



where (pf^^ and A^''^^ are the scalar and vector potentials of the nucleus at the position 
of the ith electron. These potentials, as discussed below, in general contain both PCTC 
and PVTV contributions. The second term, V^nf^\ is the PVTV contribution due to 
interaction of electron EDMs with the PCTC electric and magnetic fields {Ef^\Bf^'') 
associated with the nuclear potentials. 



Again, the "+ ■ ■ ■ " indicates additional PVTV non-electromagnetic exchanges that 
might exist between the electrons and nucleus. 

3. The interaction of the electrons with the apphed ("external") potentials {(t)f^*\ A^^^^) 
and corresponding fields (^Ef^^\Bf^^^) 



Hi$t — Kxt + Kxt ) (8) 

a (<^f - a, • J^r^) , (9) 

i=l 

(<T,-E'f^'^+ioc,.B'f^'^) . (10) 



where 



i=l 



4. The interaction of the nucleus with the external potentials 

H^^^^aj d'y [p^^\y) <t>^^^'\y) - j^^\y) ■ A^^^y)] , (11) 

where the hat over the nuclear potentials indicates that they are operators rather than 
c-number functions. Again, the nuclear charge and three-current operators include a 
variety of terms that can contribute to PCTC and PVTV interactions. A multipole 
expansion is helpful in separating the PCTC and PVTV terms. 



In what follows, we factor out electric charge from the scalar and vector potentials. 
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The spherical multipole expansion 

To divide H^^'' and H^^^^ into PCTC and PVTV interactions, we decompose the poten- 
tials, 

0W(x)=/Ap«'. (12) 

J \x-y\ 

a.A(^)(a.) = / A^^i^, (13) 
J \x-y\ 

into multipoles. One expands the photon Green's function in terms of spherical harmonics 

47r 



^ J > 



\x-y\ ^ 2 J+1 

I ^1 J>0 



>j(£)0ij(y), (14) 



where 

Y,(x) Yj{y) ^ ^f'iv) ■ (15) 

M 

The resulting expressions for the potentials are 



x-y\ 

- E ^ i ^^(^) + ^^(^)] ' (16) 

J>0 



J>1 

Here [A ® is the standard notation for the coupling two spherical tensors A and B to 
a tensor of rank J; and the "H — • " indicate the contributions from the transverse electric 
multipoles that will not contribute to either PCTC or PVTV moments due to their P and T 
transformation properties. The nuclear vector potential A*^-^^ satisfies the Coulomb gauge. 

The charge {C and C) and transverse magnetic (M and M) nuclear multipole operators 
are defined as 



= / A/y/(y)p(2/), (18a) 

= j d'y[y'Yj{y)®j{y)]^ , (18b) 

= j dhe{y-x)[{x/yf'^'-l\ y'Yf{y)p{y)^C!f^{x)-C!f^{x), (18c) 

Mfix) = J d'ye{y-x) [{x/yf-'+'-l] [y' Yj{y) ®j{y)][ 



Mf 



Cf{x 



For example, the J = 1 transverse electric multipole moment gives the nuclear anapole moment [28] . 
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Multipole Moment 
Cf. Cf{x) 
Mf, Mf{x) 
Ef, 8f{x) 



PCTC 
even J 
odd J 



PVTV 
odd J 
even J 



PVTC 



odd J 



PCTV 



even J 



Table I: Parity and time-reversal characters of the Coulomb (Cj^), magnetic (Mj^), and transverse 
electric {E^) multipole moments. 



^MUx)-MUx) 



(18d) 



In the definitions of C{x) and M.{x) - the penetration terms that account for contributions 
when the electron cloud is inside the nucleus - the subscripts "< " and ">" refer to the 
parts containing (x/y)'^'^'^^ and 1, respectively. 

The moments (diagonal matrix elements) of multipole operators in Eqs.(16-18d) have 
definite parity and time-reversal characters depending on their angular momentum, as in- 
dicated in Table I. For completeness, we have also included the moments of the transverse 
electric multipole operators, and S^{x). These are not relevant in PVTV calculations in 
which only nuclear moments are considered. A PVTV nuclear moment can arise if and only 
if some source of parity and time-reversal violation is introduced. Thus nonzero moments 
corresponding to odd charge multipoles or even magnetic multipoles are the signatures of 
PVTV interactions. 

The electron-nucleus interaction can be divided formally into PCTC and PVTV compo- 
nents: 



H 



{eM) 



int 



- OLi 



' mt — mt ' int ' 



int 



(19) 



where V^^^^ was defined in Eq. (7). In general the multipole decomposition is helpful in 
separating the first two terms. For example, in the case most often of interest in which a 
nuclear moment is being taken 



{eM) diagonal 



int 



-a J2 



An 



+ E 



J>l,odd 



.J>0,cvcn 

4:71 



2 J + 1 



Yj{xi)Q[Cj + Cj{xi 



[Yj{xi) ® a]j [Mj + Mj{xi)] 



(20) 



^7-{eAf) di agon al _ \ ^ 
''int ^ 2-^ 



i=l 



.J>l,odd 

47r 



1 



2 J + 1 
1 



Yj{xi)Q[Cj + Cj{xi)] 



J>2,even 



2 J+1 



[Yj{xi) (^a]jQ[Mj + Mj{xi)] 



(21) 
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Figure 1: The contributions to a nuclear EDM from (a) nucleon EDMs, (b) PVTV exchange 
currents, and (c) parity admixtures induced by PVTV A^A^ interactions. 



nuc 



N 



nuc 



(a) 



(b) N 




(c) A^ 




+ c. c. 



Thus one can identify the total electron- nucleus PVTV interaction V^^^^~^^''^\ 
terms where the PVTV coupling is either on the nucleus or on the electrons 



int 



int int 



the sum of 



(22) 



In connection with Eq. (21), it might be helpful at this point to comment about nu- 
clear sources of PVTV interactions, which can originate from direct PVTV electromagnetic 
couplings to the nucleon (that is, one-body PVTV currents generated by nucleon EDMs), 
couplings to various meson exchanges or to NN excitations (two-body PVTV currents), or 
PVTV interactions of any type between the nucleons. The latter would give rise to small 
PVTV admixtures in nuclear wave functions. Thus we see, working to first order in PVTV, 
that there are two classes of contributions to PVTV nuclear moments (see, e.g., Refs. [30, 31] 
and references therein): 

• The PVTV current contribution: small components of the multipole operators due to 
PVTV contributions to one- or two-body charges or three-currents - operators that 

might be denoted Cj=odd and Mj=even - would have nonzero matrix elements between 
the dominant PCTC component of the nuclear wave function. 

• The PVTV nuclear polarization contribution: large components of the multipole op- 
erators due to the ordinary PCTC currents would have nonzero matrix elements con- 
necting PCTC components in the bra state with PVTV "polarization" admixtures in 
the ket state, and vice versa. 

These contributions to the nuclear EDM are illustrated in Fig. 1. 

In the multipole expressions (Eqs. (16-21)), we have been careful to isolate terms asso- 
ciated with the finite nuclear size (and thus electron penetration) from those that would 
persist in the point-nucleus limit. Nuclear matrix elements of the operators Cj's and Mj's 
correspond to the standard static moments [29] defined by Taylor expanding the Green's 
function around = 0, a procedure that implicitly assumes |a3| |^;|. As indicated by 
Eq. (14), the complete expansion of the Green's function contains terms corresponding to 
both I a; I > \y\ as well as \y\ > \x\. Thus, to express the potentials in terms of the static 
nuclear moments, we have used 6{x — y) = 1 — 6{y — x) and collected all terms proportional 
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to 9{y — x) in the "local" multipole operators Cj(a;)'s and Aij{x).^ Physically, the static 
multipoles characterize the interaction of the electrons with a point-like nucleus, while the 
local multipoles correspond to the penetration of the electrons inside the nucleus. 

The multipole expansion also provides a useful framework for decomposing the interaction 
of the nuclear charge and current with an external potential. In the present instance we 
distinguish the terms in i?e^^ according to the transformation properties of the various 
multipole components: 

i^i^^ = + Klf ^ , (23) 

where V^^'' is even under P and T while V^^^ has the opposite transformation properties. 
For example, V^^^'' contains the interaction of the nuclear magnetic dipole moment with 
constant applied magnetic field, while V^^^ contains the interaction of the nuclear EDM 
with constant applied electric field, E^^*^ . Substituting the potential for the latter 

<t>^^^'\y)^-aEt'^ -y (24) 

into Eq. (11) leads to 

V^S^ = -ad^.Et^\ (25) 
where dj^ is the nuclear EDM given in terms of C^^ as 

dx={^-^Y Y.^^^M. (26) 

where Bm is a spherical unit vector. 

The unperturbed and perturbed Hamiltonian 

We are now in a position to identify the components of unperturbed and perturbed atomic 
Hamiltonian. The unperturbed Hamiltonian is 

z 

Ho^Yl (/^^ + a. • Pi) + l^if + + HZ'' . (27) 

i=l 

This Hamiltonian contains the following interactions: 

1. The free electron Hamiltonian, Pi me + cti-pf The electrons are treated as relativistic, 
point-like particles. 

2. The PCTC e-e interaction, V^^^^: As electrons are treated relativistically, both 
Coulomb and Breit interactions (the latter is not electrostatic) are included. The 
scalar and vector potentials in Eq. (4) acting on the ith electron in the Coulomb gauge 
are 

# = E^. Af = ^^U + .,^), (28) 

^vllGrG X i^j — I nJCj^j I — I CC^ CC j I • 
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A "local" multipole means it can only interact with another field when they overlap. The classical example 
is the extra 5 function term introduced in the Cartesian multipole expansion [29] . 
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3. The PCTC e-M interaction, V^^^'': In atomic physics, the nucleus is typically consid- 
ered as a stable, external c-number source of electromagnetic fields acting on the elec- 
trons. Here, however, we take the PCTC nuclear potentials and A^^^ [Eqs.(12- 
17)] to be dynamical quantities and expand them in terms of nuclear multipole moment 
operators, as described above. We will later resort to the c-number source approxima- 
tion when considering practical atomic calculations. 

4. The PCTC internal nuclear Hamiltonian, H^^^: Again, in typical atomic computa- 
tions, Hq is diagonalized using product wave functions that separately diagonalize 
Hf^^ and the electronic operators in Eq. (2). We will not resort to this approximate 
diagonalization until Section III, so we must include explicitly. 

The perturbed Hamiltonian contains PCTC interactions between the electrons and nucleus 
with the external field (no tilde) and PVTV interactions within the atom or with the external 
field (terms with tildes): 



T7(e) I t7(6A/') I iri.ee) 
^ext "T •'int "T ^int 



+ 



•^ext "T •'int "T ^int 



(29) 



where we have collected together the terms proportional to the electron EDM as well as the 
terms involving PVTV nuclear couphngs. The PVTV internal nuclear potential, V-^^, is 
included in the latter as one of the perturbations. Its presence will induce mixing among 
nuclear states of opposite parity. As discussed above, this mixing leads to non-vanishing 
matrix elements of the PVTV nuclear moment operators involving the ordinary nuclear 
charge and current operators. If V:^'^ carries a momentum dependence, the nuclear conti- 
nuity equation also requires the presence of exchange charge and three-current operators, 

p{y) and j(2/), that have the opposite parity and time- reversal transformation properties 

compared to the ordinary charge and current operators. Insertion of p{y) and j{y) into the 
PVTV nuclear multipole operators yields operators that have non- vanishing matrix elements 
between same-parity states. 



B. Schiff screening and its corrections 

The observabihty of the EDM of a neutral composite system is severely hampered by 
the screening effect, the cancellation between direct interactions of the external field with 
electron or nuclear EDMs and the induced terms involving the internal fields induced by 
the applied field. The terms we will discuss in the next two subsections are illustrated if 
Fig. 2. Compared to (c), (d) not only has a suppression due to the higher nuclear excitation 
energy: A£Jatom/A£Jnuc ~ 10~^, but also the e-M electromagnetic interaction here can only 
go through higher nuclear multipoles like Cq{x) and C2 or Mi, which, in combination with 
the charge dipole transition matrix element, results in an additional finite-size or hyperfine 
suppression factor, ~ 10~^ or 10~^, respectively. Therefore, the panel (d) can be safely 
ignored in the subsequent discussion. 

As Schiff pointed out in deriving his theorem [11], there exist three types of effects that 
contradict the assumptions of the Schiff theorem and, thus, make the shielding incomplete: 
(a) the constituent particles are relativistic; (b) the constituents have finite size; or (c) there 
exist non-electrostatic interactions between the constituents. All three effects can be present 
in atomic systems: (a) the atomic electrons may be relativistic, especially for heavy atoms; 
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Figure 2: Direct interactions between the external field and the electron (a) and nuclear (b) EDMs 
and the induced terms involving electronic (c) and nuclear (d) polarizations (by the external field) 
that result in screening cancellations. In the remainder of this paper, diagram (d) is ignored. 

(a) e nuc e nuc (b) 



(|AAA/^ 



e 
e 



nuc 
nuc 



e nuc 
e nuc 



+ c. c. 



Pt 



+ c. c. 



(c) 



nuc 



nuc 



(d) 



(b) the atomic nucleus has a finite spatial extent; and (c) the e-e and e-A/" electromagnetic 
interactions contain current- current components that are not electrostatic. In this section, 
all of these factors will be gathered together, and a more unified and consistent derivation 
of corrections to Schiff screening effect will be presented. 

In the literature one can find several ways of demonstrating that EDMs cause no first- 
order energy shift, such as representing an EDM as an infinitesimal displaced charge with the 
help of translational operators as in Schiff 's original paper [11], or representing the shielding 
as a hypervirial theorem with the help of the Hellmann-Feynman theorem as in Ref. [18]. In 
this paper, we follow the approach taken in Refs. [13, 14] and rely on perturbation theory. 
We concentrate first on corrections to the screening of de and follow it with a treatment of 
the more involved case of dj^. 



Shielding of Electron EDM 



To evaluate the consequences of shielding for an electron EDM, we consider the first- and 
second-order energy shifts that depend on both de and E'^^^^'^: 



^^(1) 



(g-S-lK2 |g-s-) 



(30) 
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n ° 

where |g.s.) and |n) denote the unperturbed atomic (electronic plus nuclear) ground and 
excited states, respectively: 

Ho |g.s.) = Eg,,, |g.s) , Ho \n) = \n) . (32) 

The first-order energy shift AE^^)^ arises from the interaction of the external field with the 

electron EDMs, while Ai^^j contains the effects due to the PCTC excitation of the atomic 
cloud by the external field and its de-excitation by the PVTV e-e and e-Af interactions (and 
vice- versa). 

Following Schiff, it is useful to express the PVTV interaction V^i'^ + V^iit^^ in terms of 
an appropriate commutator with Hq minus correction terms. One has ^ 

z 

^if + = - 5] [4 /3 • V, , //o] + 2 4 /3 i 75 [p1 + « {Af^ - V2 A^^) ■ p,] . (33) 

i=l 

In the non-relativistic limit, the first term on the RHS of Eq. (33) contains the commutator 
of Hq with Schiff 's displacement operator, rfg cr ■ V. Letting Hq act on either side of the 
relativistic form of this operator in Eq. (31) leads to the energy difference Eg ,. — En that 
cancels the corresponding energy denominator. One may then carry out the closure sum on 
the states leading to 

z 

A£;gj = J2 (§•«•! [deP<T,-v,,v}:i\ |g.s.) 

i=l 

+^^deYl J, I I E 75 + a ( ) - 1/2 Af ) • \n) 



E - E . 



x(n|KS|g.s.) + c.c. }. (34) 



After performing some Dirac algebra on the first term on the RHS of Eq. (34) and adding 
AE'gj and ^Ej^^}^ we obtain 

A£;g + A^} = AEl% + AE[S^ , (35) 

where 

z 

A£;{5) = (g.s. I 5^ 2 a 4 /3 ^ 75 Af"*^ • P^ |g.s.) , (36) 

i=l 



^ A similar commutation relation was also considered by Lindroth et al. in Ref. [16] and was referred to as 
"stratagem 11". It is in contrast to the earlier work by Sandars [13] ("stratagem I"), where a commutation 
relation without the (3 matrix was employed. The difference between the present work and Ref. [16] is 
that neither nuclear nor external vector potentials were included in Ref. [16]. 
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n ° ^ 1=1 

xHK2|g.s.) + C.C.I. (37) 



The effective residual first-order contribution i^E^^^^ arises from a partial cancellation be- 
tween Ai^l^j and the commutator term of ^E^!^s^ in Eq. (34). Both Ai?|^,^j and ^E^^^-^ involve 
matrix elements of 75 that connects large (upper) and small (lower) components of Dirac 
wave functions. Thus, for non-relativistic (NR) electrons (with vanishing lower components), 
Ai?|^,^^ = Ai?|2''') ~ 0) implying complete shielding of electron EDMs. 

It should be noted that the procedure for using a closure sum to evaluate the effects of 
shielding is not unique. The first discussion of relativistic effects by Sandars [13] relies on 
expressing (3 matrix in the EDM interaction as 1 + {(3 —1), and using closure just to remove 
the "1" term while leaving the (/? — 1) term intact. This approach leads to a less complete 
cancellation of AE^^j. However, as the use of {(3 — 1) also requires nonvanishing small 
components, it has the same NR limit as the 75 formalism above. In a later publication [16] , 
Lindroth et al. argued that although the 75 formalism has the advantage of being purely one- 
body when no Breit interaction is present, it is not formally of order Z"^ in contrast to the 
{13 — 1) formalism. This might introduce complications when complex atomic many-body 
calculations intend to reach this level of accuracy [16]. However, no detailed calculation 
further substantiates this claim. In this paper, the 75 formalism is chosen because it is 
simpler. 



Shielding of the Nuclear EDM 

To evaluate the effects of shielding of a nuclear EDM, we follow a procedure similar to 
that used in the case of the electron and consider energy shifts that are linear in dj^ and 
the external field to the second order in perturbation theory. However, we depart from the 
previous notation somewhat in that we make the polarization effect of V^nt"^ explicit in the 
nuclear wave functions. This yields 

Ai?f^) = (g.s.^l yif ) [g.s.^) , (38) 

For example, if the atomic wave function has the direct-product form \n) = |njv) ® [ng), 
then 

I^A^) = E F _F -^^"'l "^^^ I^^^AT ® l^)e = \^)n ® \^)e ' (^0) 



It is understood that PVTV effects are to be evaluated only in first order. Thus the meaning 
of the notation in Eqs. (38) and (39) is: 

• Contributions involving PVTV charges or three-currents in v}^^ or V^^^^ are to be 
evaluated with unperturbed wave functions. 
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• Contributions involving correspond to matrix elements with PCTC charges or 

currents in V^^^"* or t^^t'^'' • Furthermore, V^^'^ is to be treated in first order, modifying 
the bra or the ket state, but not both. 

Below we suppress the explicit subscript J\f on bra and ket states involving operators coupling 
to the nucleus, but implicitly V^^'^ has been absorbed into nuclear wave functions and its 
effects will be retained to first order in the PVTV. 

To proceed, we first write dj^ in a spherical basis using Eq. (26) and observe that 



X 



4:7l\ 1 



(41) 



Letting 



^ 1 

AHo = Ho- H^-^{Co only) = + Za V - 

i=l 



(42) 



be the unperturbed Hamiltonian without the static electron-nucleus PCTC monopole 
Coulomb interaction, we have that 



— a^. 



(Ci) 



—a 



a 



i=l 



Z 



X7 



1=1 



(43) 



Substituting Eq. (43) into A£'|^'' and letting Hq act on either side of the displacement 
operator leads to the energy difference Eg s^ — that cancels the energy denominator before 
and allows us to carry out a closure sum for this term. Doing so and letting 



int 



int 



+ a I 



(44) 



leads to 



E - E 



{(g-s-l 



1 



(g-s-l 



.s.) + c.c.j 



(e) 
ext 



1=1 



+E 



E - E 



int 



E 



^1 ■ V, , Affo' 



|g.s.) + c.c. 



(45) 



Now observe that 

z 



1 



ext 



j=l 



^3 "i ■ ^3 
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i=l 



d^f-^i, oci- A. 



(ext) 



(46) 



where we have used the electrostatic stabihty condition Z e'^'^*^ — ^y^^-^i-i _ q_ Adding 
the first- and second-order energy shifts leads to 

z 



-\{ext) 



g.s. -"^n 



i=l 



|g.s.) c.c. 



(47) 



Equation (47) demonstrates the cancellation of all terms proportional to dj^ ■ E^^^^ due 
to the screening effect and allows us to identify systematically all corrections to this screen- 
ing. We elaborate on these corrections, including the correction associated with the Schiff 
moment, below. 



Shielding Corrections 



Prom Eqs. (36, 37) and (47), one may determine all the corrections to Schiff screening 
that are linear in de or dj^ and that occur to the second order in perturbation theory. For 
purposes of future discussion, it is useful to classify the corrections as follows: 

Ground state matrix elements. The following corrections arise solely from the interaction 
with the external vector potential Af^^ : 



(1) 



1 ^ext + ^ext 



•s.), 



where 



62 = E2«rfe/3^75Af^*^p,, 

i=l 



(ext) 



(48) 

(49) 
(50) 



1=1 



Internal excitations. The remaining corrections involve excitation of atomic states by 



(e) 
cxt 



and de-excitation by operators proportional to d^, dj^, or higher PVTV nuclear moments. 
The corresponding energy shift is 



(2) 



E 



Z 9 



E - E 



(g.s.| kS \n) {n\ 62 + E E lg-^-) + ' (51) 



=1 fe=i 



where 



1=1 



(52) 
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and where, with all electronic subscripts "i" suppressed, 

. even f I — 7 ' — 

J>2 I, 

~\IjjT3 ^'+'^'''> ® ["^-^ ® Cj>{x)]j+i + {xYj{x) V)Q[dj^ , Cj{x) 

4 TT /I 4 TT Q! < — > . 

— ^ W - Y^{x) © [cZa/- Co>{x)l + {xYo{x) V)Q[dj^ , Co{x)] , 

W^Mj<{x)]j+i 

W rrnuc 

^ • V , -Hint 



of) . 



odd , 
odd . 
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even . ^ 
even . ^ 

^ = - E ^ jj-i [^•^(^) ® •^■^(^) ■ (61) 

J>2 

We refer to Appendix A for more details about the derivation of Eqs. (53-61) from Eq. (47) 
and note important conventions here: 

1. The nuclear composite operators involving dj^ and any PCTC nuclear multipole op- 
erator Xj should be realized as 

[dM X,]f = [dM ® '"^^ = 5^ {l\,jm\JM) {< , Xp}/2 . (62) 

A,m 

As two nuclear density operators not necessarily commute, the part which docs not 
commute is treated explicitly in the terms in Eqs. (53-61) involving the commutator 
of djiyf with other nuclear operators. 

2. The special gradient operator V only acts on the electronic bra and ket states as 

e ( \Oix) V I), ^ e ( , V}| )J2 = (e ( \0{x) V |), - e ( \V 0{x)\ )J/2 , (63) 

where 0{x) denotes a generic operator that depends on electronic degrees of freedom 
and where in the second equality we have performed an integration by parts. 

3. The operator "0" denotes a double scalar product as 

\,M 



We make several observations about the list of operators O^^^ ■ 

(i) The operators 0^1 are generated by the commutator of the displacement operator 
dj^f ■ V with the non-Co part of the unperturbed Hamiltonian, AHq. 

(ii) Physically, operators O^/^^ correspond to "displacement" of the static PCTC Coulomb 
and magnetic electron-nucleus interactions due to the rearrangement of the atomic 
electrons needed to maintain electrostatic equilibrium in the presence of the external 

field. The operators O^'^ describe the corresponding effects of displacing the PCTC 
penetration (non-static) e-J\f multipole interactions. 

(iii) Operator describes the response of the internal nuclear degrees of freedom to the 
external field, again as needed to maintain electrostatic equilibrium. 

(iv) Operators Oj^l characterize the effects of the "local" EDM (i.e., Ci), PVTV magnetic 
and higher (J > 3) PVTV charge multipole interactions, including both the static 

multipole interactions, , and penetration terms, ■ 
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(v) The leading, non- magnetic correction to Schiff screening due to the finite spatial extent 
and internal structure of the nucleus is the Schiff moment. It therefore arises from the 
penetration or local multipoles Cj>,<(a;) that appear with a Yi{x); a gradient acting 
on the electronic coordinates coupled to a Yq{x) or 12(2;); and the internal nuclear 
interaction in conjunction with a gradient acting on the electronic coordinates. Such 

terms appear in O'f^^ 7 . In the following section, we assemble these terms into the 
Schiff moment operator and compare with the form of the operator previously used in 
the literature. 



III. THE SCHIFF MOMENT OPERATOR 



To arrive at the Schiff moment operator, we collect all terms in O^^^^j that are propor- 
tional to Yi{x) or that contain terms with a gradient acting on the atomic coordinate and 
transform like a spherical vector: 



^1,2,5,7 



Ana 1 



a/2 

I Co>(a;) ^C2<(a;) 



ZY^{x)Q-Ci{x) 



+ 



xYo{x) V)Q[d^f, Co{x)] - {xY2{x) V)0 



dAf,-^{C2+C2{x)) 
O X 



int 



/S atomic 

= '^Schiflf 



(65) 



Note that the terms containing the gradient V acting on the electronic coordinate will 
connect electronic states differing by one unit of orbital angular momentum, as can be seen 
by performing an integration by parts and letting V act on the initial and final electronic 
wave functions. Consequently, we include these operators in 

The operator Og*°"g'^ generates the leading contributions to an atomic EDM that depend 
only on the nuclear charge density operator and not on the current density operator. The 
form given in Eq. (65) is general and could, in principle, be used to compute the atomic 
EDM generated by P- and T-odd interactions in the nucleus. Prom Eq. (51), we observe 
that 0^1°^'^ will induce mixing of opposite parity states into the atomic ground state, thereby 
allowing for a shift in the atomic energy when a constant external electric field is applied, 
i.e., Vgxt • practice, the typical atomic ground state is taken by be a direct product of the 
nuclear and electronic ground states, as in Eq. (1): 



|g.S.) = |g.S.)^(8 |g.s. 



(66) 



It is particularly interesting to consider systems in which the electronic ground state is a 
S-state (L = 0) that has a comparatively larger penetration probability due the lack of the 
centrifugal barrier. Since O^^^^ transforms as a rank-one tensor in the space of electronic 
coordinates, it will mix P states into the electronic ground state (relevant for many, but 
not all atomic EDM experiments). Thus, it is useful to derive a form for this operator 
applicable to this situation. In doing so, we can derive a nuclear Schiff moment operator 
that corresponds to - but differs in form from - the Schiff moment operator used elsewhere 
in the literature. To that end, we begin by considering the first three terms in Ogchig"^ and 
compute their S-P electronic matrix elements. 
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The electronic wave functions for the S and P states are: 

M^) = {^\S, m) = us{x) Fo(£) xT/2 , (67) 

Mx) = {x\P,\m') = up{x) Y,\x) xii2 , (68) 

where A is the magnetic projection of the P state (A = for the S state), and xY/2 ^^'^ 
spin wave function. We will henceforth neglect the spin degrees of freedom as they are 
not relevant to the derivation of the effective nuclear operator - m should be a conserved 
quantum number for this case. In order to arrive at the effective nuclear operator, we 
consider the following polynomial expansions of electronic radial wave functions near the 
origin: 

us{x)^J2^'^^''^ (69) 

fc>0 

up{x)^J2hx''. (70) 

k>l 

We now consider the matrix element of Og*^^'^ appearing in Eq. (65): 

{n\oiliTk-^.) =M {gTs.l ® (P, x\d-g^T\s) ® \grs.)j^ , (71) 

and note that we use the convention of Eq. (40). Focusing on the first term in O^^^^'^ that 
contains Co^{x) we obtain an electronic matrix element 

{p,Moii^''\sh.tter^ = -^^ I d?xrp{x)^^Mx)Y,Yr{x)d7Co>{x) 

m 

= r(^^«p(a^)*«5(a;) / dn,Y,^*{x)Y{-{x)Yo{x) 

xd]}* j d'yeiy-x)p{y)Yoiy). (72) 



The angular dQ^ integral gives Smx/V4:n. To evaluate the radial integral we first write 

up{xy us{x) =^Ckx'' , (73) 



k>l 



where the Ck are given by the appropriate products of the and bk (our procedure here is 
similar to Ref. [20]). Interchanging the order of integration and evaluating the integral over 
the electronic radial coordinate gives 

{p, A|Oslhri^)ist te.. = -^11 ^ < / ^'yp(y^ E ttt ■ ^^^^ 

Evaluating the second and third terms in the matrix element of Oschiff^ ^ similar way 
leads to 
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{k+^ Lx* 2(/,- + l) ^v^^^iAk 

- [^^ (FPI)— ^"^^ ® j ^ • ^''^ 

Note that in this expression is the nuclear operator defined by Eq. (26). 

The "local" nuclear Schifi^ moment operator can be defined by requiring that 

-{Ana) {P,X\Sl-VS^'\x)\S) ^ (P, A 101^/1^) terms 1-3 ■ (76) 
Evaluating the left side of Eq. (76) yields 

-(47ra) {P,X\S^-Vd^'\x)\S) = {47Ta) (^^^ ci 5^* . (77) 
Including for the moment only the first three terms in O^^^^^ as in Eq. (75) gives 

/k-\-l ( 

Retaining only the leading-order term, i.e., k — 1, on the RHS of Eq. (78) yields our result 
for what is conventionally referred to as the Schiff moment operator, S^: 

=^ / d'ypiy) y"[y'-i^ - ^ ® ^2(y)]? j I . (79) 

Note that nuclear matrix elements of 5^ and will involve nuclear matrix elements of 
of density-density correlations, since the last two terms on the RHS of Eqs. (78) and (79) 
contain the operator 

d^=J d?zp{z)z\ (80) 

that multiplies the operators p{y) xp' or p{y) y"^ Y^{y). 

The form of the operators in Eqs. (78) and (79) differ in two important respects from the 
operators previously used in the literature: (a) the presence of the final term involving the 
Y2{y) and (b) the present treatment of djj- as an operator rather than as a c-number. To 
illustrate the importance of this first difference, we consider the Schiff moment operator of 
the deuteron in the limit where only the nuclear polarization term (due to V^^) is retained. 
That is, we ignore the effects of PVTV one-body {e.g., the neutron and proton EDMs) and 
two-body charges and three-currents. Polarization terms (pol.) contribute via the one-body 
PCTC nuclear charge operator: 

p{y) = E ^^'^(y - ' ^P = U^+ ^3) , (81) 
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where the sum is over both nucleons {k = 1,2) with coordinate r^. Since only the proton 
contributes, we label the proton coordinate r for simplicity. Substituting Eq. (81) into 
Eq. (79) leads to 

sr^..o. =Yo-y { - + ^ ® ^^(^)]?} ' (82) 

where we have performed the integrals over y and z. Using 

[r®Y,{ft^-^r\ (83) 

we obtain 

Sr^t.ron = -lryr\ (84) 

where the coefficient —1/5 is a factor of three larger than had we omitted the l2(y) term. 
Generalizing this argument to an arbitrary nucleus we obtain 

Sr' =^E'->^{'-i-^E'^' {r}-*-^lr,0Mht^Y (85) 
Focusing on the terms in Eq. (85) for which j — k and using = gives 

where the final 4/5 arises from the Y2{y) term. For heavy nuclei, the effect of this correction 
will be suppressed by the factor of 1/Z but its importance relative to the second term in 
Eq. (79) is of the same order. The ground-state expectation value of this odd-parity operator 
will be nonzero because of PVTV admixtures in the nuclear wave function. 

The presence of the terms with j ^ k point to the impact of treating dj^ as an operator 
rather than as a c-number. To see this, we write as follows: 

S'-^ J d'y J dhp{y)y'!^y'S('Hz)-^p{z) (^z' - ^ [z ^ Y,{y)]i 

(87) 

The product of charge density operators appearing in Eq. (87) implies that nuclear matrix 
elements of involve a two-body correlation. A comparison of this result with one where 
djiyf is treated as a c-number can be made by inserting a complete set of intermediate nuclear 
states 

Ar(g:^.|p(2/)p(^)|g^.)^= J] Ar(g^.|p(2/)|n)(n|p(z)|g:i.)_^ . (88) 

n 

Treating dj^ as a c-number amounts to the assumption that this sum is effectively saturated 
by retaining only \n) — |g.s.). The numerical validity of this assumption is not at all clear. 
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and certainly should be explored in numerical calculations. This task goes beyond the scope 
of the present paper, but will be explored in future work. 

Additional contributions to the nuclear Schiff moment arise from S-P matrix elements of 
the last four terms in 0^1^^^. It is useful to distinguish contributions from these terms that 
are sensitive to electronic penetration inside the nucleus from those that are not. Following 
similar arguments to those used in deriving we obtain 

{P, MO^lutlS^^rras 4-7 = " (4 TT «) (P, A| A5(i) ■ 6^'\x) V + A5(2)0(^ Y,{x) V ) \S) 

+ (P,A|A^(3)- V |5) , (89) 



where 



, (j + 2-k) . 2V2n(j-l-k) . 



ASL = - — 



d^,C2 , (91) 



(92) 



Details pertaining to the derivation of AS^^^^ are given in the Appendix B; as for A.S^2 3)i 
they can be readily read by comparing Eqs. (65) and (89). 

The operator ^S^^^ characterizes corrections to Schiff screening that arise when the 
nuclear EDM operator does not commute with local (penetration), even- J moments of the 
charge operator. The result is the presence of the commutator [p{z) , p{y)] in Eq. (90). We 
expect that contributions from this commutator will be suppressed by powers of v/c, where 
f is a typical nucleon velocity, since the leading-order (in v/c) parts of the one-body charge 
operator give a vanishing commutator. Inclusion of sub-leading one-body or two-body charge 
operators will not, in general, give a vanishing commutator since these sub-leading terms 
contain momentum and/or Pauli spin operators. Again, we defer a detailed analysis of the 
numerical importance of these effects to a future study and concentrate here on delineating 
the various contributions to the operator. 

The operators AS'^2 3) characterize long-distance screening corrections that occur when 
the nuclear EDM operator does not commute with either the quadrupole moment operator 
or the internal nuclear Hamiltonian. We expect the former to be suppressed by powers of 
v/c for similar reasons as for ^S^iy The commutator of with Hl^'^ is, perhaps, more 
subtle. In general, the nuclear Hamiltonian contains both momentum- and spin-dependent 
forces. The commutator of the momentum and spin operators will not commute with dj^ 
which generally contains both spin- and space-dependent components. As with ^S^^^ we 
defer a detailed analysis of these corrections to a future study. 



Magnetic corrections 

Additional corrections 

As with the Schiff operator terms, we collect the lowest rank terms which can also induce 



Additional corrections to Schiff screening arise from the magnetic operators in 0^\q. 
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electronic S-P transitions: 



O 



{N) 
3,4,8,9 



- lowest rank 



-[y2(^) ® a]2 



[F2(£) ® a]i 




8 



+ 



[yo(^) a]i 
[y2(^) a]2 + A^2(a:;)) 



+(x 0a]i V)0 



Aatomic ((YX\ 

= '^Schiflf-mag ■\^'^) 



The operator Og^,^'L„^g_g contains two types of terms: (a) those arising from the "displace- 
ment" of the magnetic dipole interaction, containing the commutator of dj^ with the mag- 
netic dipole operators, and (b) an explicit magnetic quadrupole interaction. Effects of the 
former type were identified in Schiff's original paper, but have generally been neglected in 
subsequent work. For atoms with nuclear spins / = 1/2, there is no magnetic quadrupole 
contribution. 

The form of Ogl^'i^g^^ simplifies when we consider the electronic matrix element, again 
under the factorization assumption of Eq. (1). To arrive at this simplification, we consider 
the symmetry properties of the relevant electronic operator 



(94) 



To proceed further, it is useful to consider the perturbed atomic states appearing in Eq. (51) 
and label them according to their angular momentum quantum numbers. We define the 
states \J'Mj = J)^ entering Eq. (51) via 



(^lKJtV^)e K0|g.s)^, 



g.s. 



En 



(95) 



where | J, J)g denotes the electronic ground state and V^^t contains the potential associated 
with constant external electric field that mixes electronic states differing in orbital angular 
momentum by one unit. The electronic component of Eq. (51) involving a generic OJ* 
operator then has the form 



(j'j|C»°|JJ)^+e (JJ|0°|J'J) 

Now using the following identity 

{a\0\h) = {T{h)\T T-^\T{a)) , 



(96) 



(97) 



where T is the anti-unitary time-reversal operator, along with the Hermiticity and time- 
reversal properties of O'^ that 



l+l-j+m 



[Yi oc]-^ , T[Yi® a]f = (-1)'-^+- [Yi a\J , (98) 
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we can show that 



{j'j\opj)^ = -e (T(jj)|0?|r(j'j))^ 

= -(-1)^+^+^ {-ly'+j+L' (J - j\o]\j' - J) 
^_(_i)^'-^+^'+^(_i)2j+j'+j+i ^jjpp'j) , 



(99) 



where the second hne follows from the phase convention for the time-reversed state (we 
use here the Condon-Shortely convention^ ), and the third line follows from the symmetry 
properties of the 3 — j symbol. Given the fact that (L + L) is odd (parity change) and 
2 (J + J) is even, we have 



{j'jmjj)^^{-iy {jjmj'j) , 



(100) 



which means the sum, Eq. (96), can only be nonvanishing for j = even. Thus, the [10,2(2;) (8) 
a]i parts in Eq. (93) yield no contribution to the induced atomic dipole moment and we 
may work with the operator 



/Aatomic 
*^Schiff-mag 




a 2 



+[Y2(x) ® a\2 Q - M2 + M2(x) 
5 




Mi+^i>(x) + 




+{x[Yi{x)^a\i V)© 



(101) 



Several comments can be made regarding this Oschiff 



atomic -i-oT-m • 
-mag ^-cim. 



1. As O^^^^l^^^ involves interactions with static magnetic moments Mi (through the 
combination with dj\f) and M2, they are not affected by the screening. Compared 
to the finite-size effects which contain a suppression factor roughly to the order of 
(nuclear size) ^/(atomic size)^ ~ fm^/ap ~ 10~^, these magnetic effects are only sup- 
pressed by the typical hyperfine factor a"^ me/mi^ ~ 10^^ (note that we have not taken 
any atomic or nuclear enhancement factor into account, just use a pure dimensional 
analysis). These could be potentially important, as already pointed out by Schiff [11] 
and others (see, e.g., [17, 19]), particularly for the open-shell atoms. 

2. The last term in Og'^Q^^^^, unlike the rest which are of quadrupole nature and need a 
nuclear spin equal or greater than 1, contributes for any nucleus with spin. In fact, it 
has been identified in the original paper by Schiff [11] (he considered the hydrogen atom 
only) and later on studied by Hinds and Sandars [32] with a more refined expression. 
The latter authors find a non-negligible contribution from this term, about 20%, to 



Another popular phase convention is the one of Biedenharn-Rose, which introduces an extra "i^" factor 
for . But the result is independent of phase convention. 
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the EDM of TIF molecule, assuming the proton EDM is the only CPV source. We 
emphasize that this term is a result of taking dj^ and the magnetic dipole operators 
Ml, All as operators: had one started with dj^ as a c- number, this term would not 
have existed at all. 

IV. PRACTICAL CONSIDERATIONS 

In the foregoing analysis, we have laid out the structure of the atomic operators that 
describe the corrections to Schiff screening and that characterize the leading contributions 
to an atomic EDM arising from the relativistic motion of the electrons, the finite size and 
internal structure of the nucleus, and magnetic interactions between the atomic electrons 
and nucleus. In doing so, we have attempted to remain as general as possible without 
making specific reference to the atomic states. In this context, the leading contribution 
associated with finite nuclear size and internal structure is given by O^^^'^ in Eq. (65). 
Given a complete basis of atomic states including the effects of P- and T-odd admixtures 
into the nuclear states of definite parity - one could in principle use Ogchiff'^ compute the 
energy shift of an atom in an external electric field via Eq. (51). 

As a practical matter, it has been the convention to specify the computation of finite- 
size corrections by considering the effect of O'^^^^ on mixing between electronic S- and P- 
electronic states, making the factorization assumption of Eq. (1) and neglecting nuclear po- 
larization corrections. In order to compare the implications of our formulation with previous 
analyses, we have derived an effective nuclear Schiff moment operator from O^^^"^ under 
these assumptions and break down the matrix element as 

(P, \\dt^i^\S) =(-47ra) (P, A| [Sl ■V5^'\x) + ASfi) • 5^^\x) V 

+ ASf^)©^ Y2{x) V]\S) + (P, A| A5(3) • V |5) , (102) 

where S\ and AS'^^_3^ are given by Eqs. (78) and (90-92) respectively. We emphasize that 
different effective operators will apply for other atomic transitions. The derivation of the 
relevant operators starting from O^^^ will follow the similar logic as in our derivation here 
for {P^y^\0'f^^\S) . One can introduce a further level of approximation by specifying the 
sums in these equations to the leading terms, corresponding to retaining only the leading x- 
dependence of the electronic wave function near the origin. Quantifying the error introduce 
by making either the factorization approximation or retaining only the leading electronic 
x-dependence is an important task for future nuclear and atomic structure computations. 

Having introduced these two approximations, we have obtained a form for (P, A | Og^^jg"^ | S*) 
that differs in several respects from the nuclear Schiff moment operator used previously 
in the literature. These differences arise primarily because we have reformulated Schiff's 
theorem entirely at the operator level, have retained the full nuclear operator-dependence of 
the finite-size correction, and proceeded consistently within the framework of the spherical 
multipole expansion. In brief, the resulting differences with previous forms of the operator 
obtained in our approach are: 

(i) The presence of the Y2{x) term and the product of nuclear density operators in Eqs. (78, 
79); 
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(ii) The presence of the commutator of miclcar density operators in Eqs. (90, 91) that 
vanishes at leading order in v/c but will not in general vanish at higher orders due to 
the presence of spin- and momentum-dependences in the nuclear charge operators; 

(iii) The presence of the internal nuclear Hamiltonian in Eq. (92) that results from in- 
cluding the internal nuclear degrees of freedom as dynamical quantities in the atomic 
Hamiltonian. 

We have illustrated the potential impact of including some of these new ingredients by 
studying the effect of the Y2{x) term in in the deuteron, where it increases the magnitude 
of the Schiff moment by a factor of three at the one-body operator level, and in the diagonal 
part of the density-density operator matrix element for heavy nuclei, where its impact is 
1/Z suppressed. We have not, however, quantified the effect of non-diagonal terms and the 
corresponding nuclear excitations, those that arise at higher-order mv/c, or those associated 
with the internal nuclear Hamiltonian. The result of future studies that quantify these 
contributions will determine the degree to which previous computations of the Schiff moment 
adequately characterize the leading nuclear correction to Schiff screening. 

Additional corrections arise from the effect of magnetic interactions between the atomic 
electrons and nucleus. In principle, these corrections could be important for the theoretical 
interpretation of other atomic EDMs. To that end, we have worked out the leading form 
of the magnetic operator O'^^^t^a^g-, again in a way that makes no reference to the atomic 
states and later specifying to the situation of a simple direct product of nuclear and electronic 
states. Analyzing the quantitative impact of this operator will also be an interesting endeavor 
for future nuclear and atomic structure studies. 
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Appendix A: SYMMETRIZATION OF [A{y) B{x) , C(y) D{x)] 

When evaluating Eq. (47) for the residual internal PVTV e-J\f interaction that evades 
the Schiff screening for dj^, the commutator • V , Aifo], can be generically expressed as 
[A(y)B(x) , C{y) D(x)] with 0{y) and 0{x) denoting operators acting upon nuclear and 
electronic Hilbert spaces, respectively. In order to obtain a manifestly Hermitian result, we 
employ the following identity: 



where we have used the fact that operators depending on x and y commute: [A{y) , B{x) ] = 
0, etc. As an illustrative example, we go through the terms in AHq which involve static 
charge multipoles Cj's with J > 2. For this case, one takes 



[A{y) B{x) , C{y) D{x)] =1/2 {[B{x) , D{x)] {A{y) , C{y)} 




(Al) 
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C{y) = Cf , D{x) = Yf*{x) , (A2) 
and uses the gradient formula 

where Yjf^ is the vector spherical harmonics, then the commutator can be reduced as 



C^AT • , Cj Yj{x) 



+ Yj(x) V 0|dAf , Cj] , (A4) 



with the short-hand notations [...]®^™, V, and 0, being defined in Eqs. (62-64). This result 
leads to the operator, i.e., Eq. (53). 

Appendix B: DERIVATION OF AS^i) 

In this appendix, we use the term involving Co{x) in Eq. (65), i.e., the 4th term, to 
illustrate the general procedure for obtaining A<S'(i). First, the electronic matrix element 

{P,MOf^''\S),,,^,^{P,X\^Yoix) [dj,,Coix)] ■ V]\S) , (Bl) 

Z/ X 

can be easily evaluated with 

/I X 
d'x - Yo{x) e{y -x) (--!)- {rp,x - ^5 V^p,a) 
X y A 

SO that 

(^,A|6sral'5)term4 

Since this specific matrix element has been defined as (see Eq. (89)) 

(P, \\df^^\S)^,^^ ^ -(4vr«) (P, A| A5(i) 5^'\x)V\S) (B4) 

= -(47ra) \haoSl^, (B5) 



2 v3 Z ^ , - „ 
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then the contribution of the Co term to S^^^ becomes 

j>l,k>0 V 1 0/ J J 



(j + A: + 2)(j + fc + l) ■ ^ ^ 



Similar logic applies to the term involving C2{x). 
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